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1.  Introduction 

We  have  now  eompleted  a  Phase  II  development  of  an  atmospherie  deeision  aid  (ADA) 
foreeasting  methodology  for  military,  eivilian,  and  eommereial  aireraft  for  whieh  signifieant 
wave  and  turbulenee  aetivity  may  pose  an  operational  or  funetional  risk.  The  speeifie  goal  for 
MDA  purposes  was  to  ereate  a  foreeasting  methodology  for  turbulenee  aetivity  at  the  expeeted 
High  Altitude  Airship  (HAA)  flight  altitude  of  22  km  that  speeifioally  addresses  all  of  the 
dominant  sourees  of  turbulenee  in  the  elear  atmosphere.  These  inelude  turbulenee  aeeompanying 
gravity  waves  arising  from  deep  eonveetion,  those  arising  due  to  air  flow  over  signifieant  terrain, 
and  turbulenee  due  to  Kelvin-Helmholtz  instability  (KHI)  aeeompanying  wind  shears  arising 
from  jet  stream  flows  and  inertia-gravity  wave  motions  in  the  lower  stratosphere. 

Given  the  expeeted  deployment  of  HAAs  over  speeifie  sites,  our  implementation  for 
HAA  purposes  provides  foreeast  options  extending  up  to  72  hours  within  a  user-seleeted  domain 
either  500  or  1000  km  aeross  eentered  at  any  point  on  Earth.  Foreeasts  are  based  on  National 
Center  for  Environmental  Predietion  (NCEP)  initial  eonditions  and  NCEP  foreeasts  at  6-hour 
intervals  thereafter.  As  with  any  foreeasting  methodology,  however,  we  have  the  greatest 
eonfidenee  in  the  foreeasts  based  on  direet  measurements  at  the  earliest  foreeast  times. 

The  motivations  for  our  approaeh  and  the  speeifie  methodologies  employed  for  eaeh 
eomponent  of  the  turbulenee  foreeast  were  deseribed  in  previous  reports.  This  final  report 
summarizes  our  development  effort  and  the  resulting  turbulenee  foreeast  produet  that 
aeeompanies  it  and  whieh  is  provided  as  a  suite  of  operational  seripts  and  exeeutable  files.  The 
operational  and  support  environment  for  the  foreeasting  methodology  are  deseribed  in  Seetion  2. 
Examples  of  the  foreeast  output  for  various  sourees  and  environments  are  provided  in  Seetion  3. 
Seetion  4  briefly  deseribes  the  required  operational  environment.  A  separate  HAA/TURBO 
Users  Manual  also  aeeompanies  this  report  and  the  seripts  and  exeeutable  foreeasting  files. 

2.  HAA/TURBO  Operational  and  Support  Environment 

The  eomputational  and  data  aeeess  requirements  neeessary  to  perform  turbulenee 
foreeasts  with  our  HAA/TURBO  software  are  deseribed  in  the  User  Manual  aeeompanying  the 
exeeutable  eodes.  The  nominal  data  feed  requires  input  from  a  NOAAport  downlink  of 
initialization  data  or  other  suitable  data  souree  providing  input  data  in  the  same  formats. 
Additional  details  are  provided  in  the  User  Manual. 

3.  HAA/TURBO  Methodologies,  Turbulence  Intensities,  and  Example  Forecasts 

a,  HAA/TURBO  forecast  methodologies 

Four  specific  methodologies  underpin  our  HAA/TURBO  turbulence  forecasting 
procedure.  As  described  in  our  previous  reports,  these  include 

1.  a  Fourier-Eaplace  (EE)  method  developed  by  Vadas  and  Fritts  (2001,  2002,  2004, 
2008)  for  the  description  of  the  gravity  waves  arising  from  deep  convection  with  fast 
embedded  convective  plumes  (Fane  and  Sharman,  2006)  and  their  propagation  to 
higher  altitudes  in  varying  environments; 

2.  a  Fourier-ray  (FR)  method  describing  transient  and  steady  gravity  wave  responses  to 
complex  3D  terrain  and  their  propagation  to  high  altitudes  (Broutman  et  ah,  2003, 
2006,  2008;  Eckermann  et  ah,  2006); 


2 


3.  characterization  of  the  statistical  potential  for  KHI  due  to  wind  shears  and  based  on 
the  dependence  of  such  instabilities  on  the  mean  and  low-frequency  wave 
environment  at  HAA  altitudes;  and 

4.  characterization  of  turbulence  intensities  for  both  gravity  wave  breaking  and  KHI 
based  on  direct  numerical  simulations  (DNS)  of  these  dynamics  that  specifically 
resolve  the  large-scale  instability  dynamics,  the  inertial  range  of  turbulence,  and  the 
turbulence  variability  throughout  the  domain  and  the  event  evolution  (Werne  and 
Fritts,  1998,  1999a,  b,  2000,  2001;  Werne  et  ah,  2005;  Fritts  et  ah,  2008a,  b). 

These  individual  methods  vary  greatly,  but  they  have  each  enabled  development  of  a 
quantitative  method  for  forecasting  turbulence  intensities  that  allow  a  merged  forecasting  scheme 
accounting  for  all  important  turbulence  sources  based  on  NCEP  initial  conditions  and/or 
forecasts.  These  modules  are  linked  in  order  to  provide  a  combined  forecast  for  mean  turbulence 
intensities  and  their  expected  hazards  for  HAA  operations  as  shown  in  Figure  1.  For  each 
module,  we  have  employed  the  best  and  most  quantitative  data  available  for  assessing  the 
individual  contributions  to  the  total  expected  mean  turbulence  intensity  at  50-km  resolution 
across  the  forecasting  domain.  Examples  of  the  turbulence  structures  for  KHI  and  GW  breaking 
from  which  we  determine  turbulence  statistics  are  shown  in  Eigure  2. 


Figure  1,  Schematic  of  the  merged  forecasting  methodology  employing  modular  assessments  of 
mean  turbulence  due  to  each  of  the  three  primary  sources  of  clear  air  turbulence  at  the  HAA 
flight  altitude  of  22  km;  mountain  waves  (MWs),  gravity  waves  (GWs)  generated  by  deep 
convection,  and  KHI  occurring  due  to  strong  wind  shears.  The  modular  approach  both  accounts 
for  all  of  the  major  turbulence  sources  and  allows  for  forecast  updates  as  new  data  are  available 
and  forecast  improvements  in  eaeh  module  as  initial  data  or  forecast  capabilities  improve. 
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Figure  2,  Examples  of  turbulence  due  to  KHI  (upper  left),  mountain  wave  or  convective  gravity 
wave  breaking  (right,  times  are  in  buoyancy  periods  and  turbulence  initiation  occurred  at  t  ~  10), 
and  the  distributions  of  logic  e  for  the  GW  breaking  event  (lower  left).  Characteristic  vertical 
depths  of  these  turbulent  regions  are  ~100  m  to  1  km  for  KHI  and  --1  to  10  km  for  GW  breaking. 
The  distributions  of  turbulence  intensities  within  each  type  of  event  indicate  a  wide  range  of 
values,  with  the  statistics  dependent  on  the  event  scale  and  environmental  factors. 


Turbulence  statistics  for  both  mountain  waves  (MWs)  and  the  convective  gravity  waves 
(GWs)  are  based  on  high-resolution  DNS  of  these  dynamics.  These  statistics  (see  the  lower  left 
panel  of  Figure  2)  indicate  a  broad  range  of  turbulence  intensities,  with  mean  values  at  each  time 
far  above  minimum  values  and  also  well  below  maximum  local  values.  Following  conventional 
usage  for  aircraft,  we  will  categorize  turbulence  as  “very  light”,  “light”,  “moderate”,  “severe”,  or 
“extreme”,  with  each  category  spanning  a  decade  or  more  of  turbulence  intensities.  Hence, 
“extreme”  turbulence  intensities  may  occur  in  localized  regions  when  mean  turbulence 
intensities  are  “severe”,  “severe”  turbulence  may  occur  in  localized  regions  when  mean 
turbulence  intensities  are  “moderate”,  and  so  on,  with  the  probability  of  local  values  greater  than 
10  times  the  mean  value  (i.e.,  “severe”  rather  than  “moderate”  or  “moderate”  rather  than  “light”) 
of  ~1%.  Thus,  we  have  designed  a  forecasting  procedure  for  anticipated  mean  turbulence 
intensities,  but  we  note  that  extrema  within  each  category  are  an  unavoidable  consequence  of  the 
natural  variability  within  any  turbulence  field.  We  also  note  that  the  potential  for  the  extrema 
within  each  category  to  impact  HAA  flight  operations  or  task  performance  depends  on  the  spatial 
scales  of  the  turbulence  and  that  this  increases  with  mean  turbulence  intensity. 
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b.  Mean  turbulence  intensities 

Measurements  of  stratospherie  turbulenee  over  many  years  have  led  to  a  qualitative 
understanding  of  the  mean  and  extreme  turbulenee  intensities  due  to  the  various  sources  posing 
the  greatest  risks  to  aircraft  in  the  clear  atmosphere.  These  include  aircraft  measurements  at 
altitudes  as  high  as  21  km  and  balloon  measurements  extending  to  similar  altitudes  (Lilly  and 
Kennedy,  1973;  Lilly  et  al.,  1974;  Cadet,  1977;  Cot  and  Barat,  1986).  Based  on  these 
measurements,  and  their  stratification  by  meteorological  conditions,  we  estimate  a  mean 
turbulence  intensity  due  to  background  shear  instability  having  a  mechanical  energy  dissipation 
rate  of  e  ~  10"^  m^s'^.  Trout  and  Panofsky  (1969)  related  energy  dissipation  rates  to  aircraft 
assessments  of  turbulence  intensities  in  the  categories  identified  above,  and  these  are  the  basis 
for  our  relation  to  potential  flight  hazards  displayed  in  Table  1.  We  note,  however,  based  on  both 
our  DNS  of  turbulence  sources  and  evolutions  and  more  recent  high-resolution  measurements  at 
small  scales,  that  turbulence  is  virtually  always  present  at  some  level.  Thus  we  have  departed 
from  the  terminology  of  Trout  and  Panofsky  (1969)  by  defining  categories  of  “very  light 
turbulence”  rather  than  “no  turbulence”  and  “extreme”  turbulence  higher  than  “severe”  levels. 


Turbulence 

intensities 

Very  Eight 
Turb.  (VET) 

Light  Turb. 

(LT) 

Moderate  Turb. 
(MT) 

Severe  Turb. 
(ST) 

Extreme 
Turb.  (ET) 

E  (m'^s"'*) 

<  3x10'" 

3x10-^- 3x10'" 

3x10'" -3x10-^ 

3xl0-^-3xl0-‘ 

>3x10-* 

Table  1,  Correspondence  of  mechanical  energy  dissipation  rate  to  defined  level  of  flight  hazard. 


The  turbulence  intensities  defined  as  “severe”  turbulence  contain  all  of  the  estimates 
listed  by  Trout  and  Panofsky  (1969)  and  those  provided  by  Lilly  and  Kennedy  (1973)  in  this 
category.  They  have  also  been  confirmed  and  further  quantified  through  scaling  of  our  DNS 
results  to  the  Lilly  and  Kennedy  (1973)  observations,  which  were  found  to  be  in  good  agreement 
for  the  mountain  wave  spatial  and  temporal  scales  indicated  by  these  observations.  We  also  note, 
however,  that  assessments  of  turbulence  intensity  are  likely  to  be  airframe  dependent.  As  such, 
specific  applications  to  HAA  may  require  further  quantification  and  validation  when  these 
systems  are  deployed.  Applications  for  other  military  aircraft  may  likewise  benefit  from  further 
evaluation  of  turbulence  impacts  on  aircraft  having  very  different  aerodynamic  responses  than 
the  aircraft  yielding  data  from  which  the  current  thresholds  were  derived. 


c.  Example  forecasts 

Our  methodologies  for  assessing  mean 
turbulence  intensities  arising  from  the  three 
sources  also  lead  to  statistical  distributions  of 
mean  turbulence  that  are  dependent  on  the 
scales  and  environmental  conditions  of  the 
individual  events.  An  example  of  the  mean 
turbulence  statistics  obtained  from  our  KHI 
turbulence  module  for  NCEP  balloon 
soundings  is  shown  in  Figure  3.  The  mean  is 
E  =  lO’"^  m^s'^,  but  the  distribution  spans  5 
decades  of  mean  turbulence  intensities. 


LoglO(cps) 


Figure  3,  Histogram  of  mean  turbulence  forecasts  based  on  --5  million  NCEP  balloon  soundings. 
Only  ~5%,  0.6%,  and  0.01%  of  these  values  correspond  to  “light”,  “moderate”,  and  “severe” 
turbulence,  respectively.  Turbulence  due  to  MWs  and  convective  GWs  is  significantly  stronger. 
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Examples  of  our  HAA/TURBO  forecast 
inputs  and  outputs  for  a  lOOOxlOOO-km  domain 
centered  on  Denver  for  a  case  in  which  wind  shears 
and  convection  contributed  only  weakly  (a  typical 
winter  response)  are  shown  in  Figure  4.  The  NCEP 
model  wind  and  temperature  profiles  and  the  terrain 
(upper  panel)  yield  the  MW  activity  displayed  in 
the  vertical  velocity  field  (middle  panel),  from 
which  turbulence  intensities  are  derived  and  flight 
hazards  are  assessed  (lower  panel).  As  noted  above, 
these  represent  flight  hazards  based  on  mean 
turbulence  values,  with  expectations  that  local 
extrema  within  each  forecast  zone  may  be  ~10 
times  larger  (thus  one  higher  category  of  turbulence 
severity)  with  a  ~1%  probability.  The 
HAA/TURBO  forecast  for  convective  GWs 
requires  considerably  more  time  to  access  NCEP 
convective  module  initialization  data  and  perform 
the  requisite  computations.  We  expect  that  these 
procedure  will  be  streamlined  with  time  and  that 
future  HAA/TURBO  versions  will  yield  forecasts 
more  quickly. 

Figure  4,  Examples  of  the  terrain  database  (top), 
MW  vertical  velocity  field  (middle),  and 
HAA/TURBO  forecast  of  flight  hazard  at  22  km  in 
50x50-km  tiles  (bottom)  for  a  lOOOxlOOO-km 
domain  centered  on  Denver  (with  green,  yellow, 
and  orange  tiles  corresponding  to  “VET”,  “LT”, 
and  “MT”,  respectively,  see  codes  in  Table  1). 


Topography,  m 


-400  -200  0  200  400 
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Vertical  Velocity,  m/s 


d.  Other  forecast  outputs 

Our  forecast  methodologies  also  permit  the 
computation  and  display  of  various  intermediate  or 
additional  forecast  products,  which  may  also  be 
useful  for  HAA  applications.  Possibilities  include 

1.  vertical  velocities  due  to  MWs  and 
convective  GWs  (as  in  Figure  4); 

2.  maps  of  energy  dissipation  rates,  e, 
from  which  turbulence  intensities  are 
assessed,  at  higher  resolution  than 
employed  for  turbulence  intensity 
categories; 

3.  maps  displaying  convective  cell 
locations,  intensities,  and  the  related 
turbulence;  and 
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4.  various  meteorological  fields,  possibly  including  NCEP  horizontal  winds,  frontal 
systems,  etc.,  surface  maps  showing  roads  boundaries,  infrastructure,  and/or  other 
user-defined  fields. 

Hence,  feedback  providing  preferences  or  suggestions  for  additional  fields  that  would  be 
valuable  to  HAA  forecasters  are  encouraged  and  will  be  considered  for  inclusion  in  future 
HAA/TURBO  software  updates. 

4.  HAA/TURBO  Operational  Factors 

Our  HAA/TURBO  software  has  been  configured,  compiled,  and  tested  on  a  Linux 
platform  having  1  GB  of  memory.  However,  we  recommend  a  minimum  of  2  GB  of  memory  to 
ensure  that  the  software  has  sufficient  memory  resources,  given  other  processes  that  may  also  be 
running  in  parallel.  We  have  configured  the  software  to  provide  what  we  consider  to  be  optimal 
resolution  for  the  convective  GW  computation,  and  this  requires  approximately  30  min  to 
execute  for  each  forecast  time  desired.  This  can  be  changed  to  require  significantly  less 
computational  resources,  but  at  some  cost  in  forecast  sensitivity  to  all  GW  spatial  scales.  Our 
graphical  display  package  currently  employs  IDL  software  and  routines,  but  can  likewise  be 
reconfigured  with  other  public-domain  software  if  desired.  These  and  other  operational  aspects 
of  the  HAA/TURBO  forecast  software  are  described  in  greater  detail  in  the  accompanying  User 
Manual. 
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